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1. Introduction 

We study nonlinear systems, describing impact dynam- 
ics. One of the most interesting phenomena, occurring in 
impact systems, is the so-called grazing, where impact ve- 
locity tends to zero. The stability of the impacting mo- 
tion was first comprehensively considered by Peterka [f]. A 
decade later Thompson and Ghaffari [2] conducted numeri- 
cal studies showing complexity of dynamical responses for a 
simple impacting system. Mathematically, a grazing bifur- 
cation corresponds to the case when there exists a family of 
periodic solutions, continuously depending on a parameter 
of the system and such that the normal component of the 
impact velocity vanishes at the bifurcation point. 

There is a rich literature on analysis of the dynamics 
near grazing in a general class of impact oscillator systems, 
most notably the comprehensive work of Nordmark [3-10]; 
see [11,12] for an overview. A powerful technique of analysis 
is to derive a so-called discontinuity mapping which has a 
square- root singularity, close to the impacting orbit. This 
can then be combined with an analytic Poincare map to 
give a so-called grazing normal form whose dynamics can 
be shown to be topologically equivalent to those of the un- 
derlying flow. See also the works of Whiston [13-15], Budd, 
Dux [16,17] and Chillingworth [18,19] for some geometrical 
analysis of the underlying strange attractors in impact os- 
cillators. The structure of the unstable manifold of a limit 
cycle near grazing was studied in the paper [19]. 

In this paper we have a different goal, namely to show 
that the topological features of the resulting attractors can 
be characterized in terms of the standard tools for chaotic 
dynamics. In particular, we use the following definition of 
chaos, given by Devaney [20] . 

Definition 1. Let U, V € M be domains of a smooth 
manifold. Consider a diffcomorphism T : U — > V. An in- 
finite compact hyperbolic invariant set K C U f] V is called 
chaotic if the following conditions are satisfied. 

1. The periodic points of the mapping T are dense in K. 

2. There is a point p G K, whose orbit 

O(po) = {T k ( Po ) :keZ} 

is dense in K. 

The periodic points of a chaotic invariant set must be un- 
stable. An attractor is called strange if it contains a chaotic 
invariant set. 



It was shown by Nordmark [3,6] that the periodic solu- 
tions must be unstable in a small neighborhood of grazing. 
The behavior of the corresponding Lyapunov exponents was 
studied. It was shown in many papers (experimentally [22- 
26], numerically [27-30] and analytically [6], [31-34] that 
the grazing is one of typical reasons for chaotic behavior in 
vibro- impact systems (VIS). 

Usually, the Devaney chaos is found via existence of a 
transversal homoclinic point, corresponding to a transversal 
intersection of the stable and the unstable manifold of a 
fixed point (the Smale-Birkhof theorem, [35] ) . 

An approach to study the structure of the stable and un- 
stable manifold was offered in the article [28] (see also [29]). 
The homoclinic point is obtained by calculation of a bent of 
the unstable manifold in the points of non-smoothness. This 
approach was applied for VIS in the current paper (Section 
4). 

The other types of chaos are also possible for vibro- 
impact systems. The sufficient conditions for the Li-Yorkc 
chaos in vibro-impact systems have been given in [11,30] 
and [36] , the stochastic properties of impact oscillators were 
studied in [37,38]. 

However, due to the strong nonlincarity of vibro-impact 
systems, especially the non-s.d.f. ones, the ratios of Lya- 
punov exponents may be very big for a near-grazing period- 
ic solution. This follows from the normal form of the square 
root singularity [6,11]. Hence, the corresponding strange at- 
tractor may be very sensitive to the changes of parameters 
of the systems. It may become invisible [25] i.e. it does not 
physically occur but could occur if the appropriate param- 
eters and initial data were chosen. 

Roughly speaking, we have an object (a part of an at- 
tractor) which looks like a periodic point on the computer 
screen. On the other hand the structure of invariant mani- 
folds of this object may be complex, they may look like fo- 
liations over the Cantor set. Then, the following questions 
arise. 

1 . Does the inner structure of an invisible attractor have 
an influence on the global behavior of solutions? 

2. How can the invariant manifolds, corresponding to 
subsets of an invisible attractor look like? 

This problem has been analyzed by one of the co-authors 
[25,26]. For a linear harmonic oscillator with an elastic im- 
pact the structure of the unstable manifold of the periodic 



2 



point was studied (Fig. 1). The following system is consid- 
ered: 

x = y; 

y = au 2 sin(wi) — 2£y — x — (3{x — e)H(x — e). 

The last term at the second equation corresponds to the 
soft impact (the stiffness of the delimiter (see Fig. 2) is 
much larger than one of the spring). Both the attractors 
at the Fig. 1 correspond to the values £ = 0.01 (damping), 
e = 1.26 (position of the delimiter), a = 0.7 (excitation 
amplitude), (3 — 29 (stiffness ratio), H is the Hcaviside 
function. The considered grazing bifurcation corresponds 
to lo = 0.801928. 

The left hand side of the figure shows an unstable man- 
ifold of a 3-periodic point, u> = 0.802 the right hand side 
shows one of a chaotic invariant set, u = 0.8023, which oc- 
curs for a very narrow band of the parameter w(see [25] for 
details). One can see that whatever the considered invari- 
ant set is, the corresponding unstable manifold has a visible 
complex structure. So, the structure of an unstable mani- 
fold or, more precisely, one of its numerical approximation, 
traces one of the corresponding invariant set, but is much 
more easy to be studied numerically. 




Displacement Displacement Fig. 1. 

In the current paper, we give an attempt to explain such 
phenomena and offer an approach, which gives a new type 
of chaotic dynamics. Instead of considering the dynamics of 
points, we introduce (Theorem 2) the dynamical system on 
an infinite set of disks of codimension 1. We show that this 
dynamics is transitive and periodic disks are dense. Con- 
structing these so-called " admissible disks" approximating 
the invariant manifolds, simplifies the studying of the inner 
structure of the invariant set. 

The obtained dynamics is structurally stable. Unlike the 
classical results of chaotic dynamics (like the Smale-Birkhof 
theorem on the chaotic dynamics in a neighborhood of a 
transversal homoclinic point), we do not assume that all 



the Lyapunov exponents of the bifurcating periodic/fixed 
point are non-zero in a neighborhood of the bifurcation val- 
ue. Comparing with the similar results from the theory of 
partial hyperbolicity [39-47] (see also references therein), 
we do not need any special assumptions like dynamical co- 
herence (integrability of the central foliation) or dimension 
1 of the central manifold. 

This dynamics of disks may be considered as a new mod- 
el of the near-grazing oscillations. 

The rest of the paper is organized as follows. At the Sec- 
tion 2 we introduce the Newtonian model of a vibro-impact 
system and discuss the basic properties of the correspond- 
ing solutions. At the Section 3 we define a grazing family of 
periodic solutions and analyze, when such a family exists. 
The Section 4 is the main. There, we provide some condi- 
tions, sufficient for existence of a non-hyperbolic homoclinic 
point in a neighborhood of grazing. The existence of such 
a point does not imply the classical Devaney chaos, so we 
introduce a special dynamical system on the space of the so- 
called admissible disks which is described by the symbolical 
dynamics. A simple example, illustrating the main result of 
the paper, is considered at the Section 5. At the Section 6 
we discuss how the obtained results may be applied. The 
conclusions are given at the Section 7. 

Some definitions and results are accompanied with re- 
marks where the corresponding statements are explained at 
a less formal level. 

2. Mathematical model of an impact oscillator. 

2.1. General assumptions and equations of motion of the 
free flight 

Let n e N be the number of degrees of freedom of the 
considered motion. Consider the half-space 

A = [0,+oo] x R 2n -\ 

Denote by | • | the Euclidean norm. Let col(ai, . . . , a m ) be 
the column vector consisting of the elements a\, . . . , a m . For 
any vector 

b = col(6i, ...,&„) 

we define b tan = (62, • • • , b n ). 

Suppose that the free flight motion between impacts is 
described by the following equations 

Xk fkif-, x\ : Xi, . . . , x n: x n: ^/), k 1, . . . , n. 
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or, equivalently, by the system 

Xk = Vk, Vk = fk(t,z,fi), k = l,...,n. (1) 

This system may be rewritten in the form 

z = F{t,z,n) = col(yi,/i(f,z,/x),.. .,y n ,f n (t,z,n)), 
z = col(a;i,yi, . . . ,x n ,y n ). 

We suppose that the parameter [i is scalar. The continuous 
function F : R 2 "+ 2 — > K 2 " is assumed to be T - periodic 
with respect to the variable t and C 2 smooth with respect 
to z and [i. 

2.2. Impact conditions 

Suppose that the system (1) is defined for z € A and 
if X\ = the following impact conditions takes place. Let 
fi- < < /i+ be a segment, r : [//_, fi + ] x OA — > [0, 1] be a 
C 2 smooth function. 
Condition 1. 

1. If the solution 

z(t) =co\(zi(t),...,z n (t)) = 
col(xi (t ) , yi (t ) , . . . , x n (t ) , y n (t ) ) 

of the system (1) is such that 

ar 1 (t -0)=0, i/i(t -0)<0, 

a;(to+0) =ar(t -0), 
tfi(<d + 0) = -r(A*)yi(to-0), 

2/tan(*0 + 0) = ytan(t ~ 0) 

where y a (to — 0) anrf y Q (io + 0) a^e velocities before 
and after impact respectively, a is 1 or tan. 

2. Let the solution z(t) of the system (1) be such that 
xi{to) = 0, yi{ta) = for a certain instant to. Let 

be the solution of the system 

Xk = Vk', Vk = fk(t, 0, 0, X t an,ytan,V), 

k = 2, . . . , n. 

with initial conditions C(^o) = z ta n{to), and a segment 

i=[to,h} 

be such that 

/i(t,o,o,c(t),/i)<o 



for all t\ € /. Then 

aji(i)=0, yi(t) = 

and Ztanit) = C,(t) for any t E I. 

Remark 1. So, we suppose, that xi is an impacting 
variable, the delimiter is homogeneous (the restitution co- 
efficient does not depend neither on the point x tan nor on 
the velocity y of the impact) and slippery (the tangent com- 
ponent of the velocity does not change during the impact). 

The item 1 of the Condition 1 describes the instanta- 
neous reflection, corresponding to a non-degenerate impact 
(with a positive value of the normal velocity) , and the item 
2 describes the sliding along the delimiter, corresponding to 
a zero- velocity impact. This part of the Condition 1 forbids 
the zero-velocity penetration through the delimiter which 
is a nonsense from the physical point of view. 

We study the following model of vibro-impact system: 

z = F(t, z, /i) 

while xi(t) > 0; (2) 
the Condition 1 is applied if x\{t — 0) = 0. 

Remark 2. One can keep in mind the following simple 
case of the system (2): 

x = y; y = — x — O.ly + sini, i,t/eM 

while x > and y(t + 0) = -y(t - 0) if x(t) = 0. The 
mechanical system, corresponding to this model, is shown 
at the Figure 2. Here 1 is a mass, 2 is a spring, 3 is a 
delimiter and 4 is a damping element. 




Fig. 2. 



The following statement is evident due to the implicit 
function theorem. 

Lemma 1. Let z(t) = co\{x\{t), yi(t), . . . , x n (t), y n (t)) be 
the solution of the system (2) for fi = fj,o and the ini- 
tial conditions z(to) — z° — col(x°, . . . , x„, j/„) where 
X® ^ 0. Suppose that this solution is defined on the segment 
[t_,t+] 3 to. Assume that there are exactly N zeros 

t- < t" < . . . < Tpf < t + 
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of the function Xi(t) over the segment [t-,t + ] and 

1/1(7-? - 0) ^ 0, (j = l,...,N). 

Then for any e > there exists a neighborhood U of the 
point col(z°, /Wo) such that for any fixed 

N 

teJ=[t^t + }\\J(r° k -e 7 T a k +e) 
k=i 

the mapping z(t, t\, z 1 , /x) is C 2 smooth. For these solutions 
there are exactly N impact instants 

^■(ii.zS/zi), j = l,...,N 

over the segment [t-,t+]. These instants and corresponding 
velocities 

Y j = -yiiTjitl^ 1 ,^) - O,*!,^ 1 ,^!) 

C 2 smoothly depend on t\,z^ and 

So, the instants and velocities of non-degenerate impacts 
smoothly depend on the initial data and parameters. 

3. Grazing in the Newtonian model 

Let /1* E (0, Suppose the following condition is sat- 
isfied. 

Condition 2. (Fig. 3) There exists a continuous family of 
T - periodic solutions 

<f(t, A*) = col(^(i, (j,), ip\(t, fi), . . . , <p*(t, n), ip v n (t, (j,)), 

t E R, n E [0,/z*) 0/ i/ie system (2) with the following 
properties. 

1. For any /1 > t/iere exist exactly N + 1 distinct zeros 
To(/z), . . . , 7at(m) 0/ f/ie component ipf(t,(i) over the 
period [0,T). 

I?. T/ie velocities yok(n) = —^p\{ T k{p) — 0,/i) are smc/i 

yoo(M) > /or aZZ /x > 0, yoo(0) = 0, 
/i(r (0), 0, 0, ^„„(td(0), 0), 0) = 0o > 0, 
ltofc(/*)>0, V^G[0, M *), fc = l,...,JV. 

(3) 

.Here 

Ptan(*, A*) = (<P2 (*) A*), ¥>"(*, A*), ■ • • , Vn(*i A«)> A*))' 



5. TTie instants Tk(p) and the velocities yokit 1 ) continu- 
ously depend on \i E [0,/x*). 

Remark 3. So we consider the family of periodic mo- 
tions, corresponding to A~ non-degenerate impacts over the 
period and exactly one impact, such that the normal com- 
ponent of the corresponding velocity vanishes for fi = 0. 
This family may persist for negative values of \i (the so 
called continuous grazing. Figure 3) or disappear (a so- 
called discontinuous grazing). In the rest of this section 
we write down the conditions, sufficient for both of these 
bifurcations. A simple example of the system with a graz- 
ing family of periodic solutions (Eq. (27)) is studied at the 
Section 5. 




We suppose without loss of generality that To(fi) = 0, 
Vtan (0j A 4 ) = 0. This may be obtained by the transformation 
s = t- T (fj,), ( = z - col(0,0,<p ton (0,/z)). Define 

#o = - rnin f min Ti([i),T — max rjv(/i) J . 
2 \»e[o,n*) A*e[o,ji*) / 

Fix small positive 9 < 9q and fi E (0, a«*) consider the shift 
mapping for the system (2), given by the formula 

S li ,e{z o ) = z{T-e + O,-0,z°,fjL). 

For small positive n and 6 the mapping S^g is C 2 smooth 
in a neighborhood of the point z^o = ip(—6,iJ,). Then 

lim z^q = 0. 

The existence of a grazing family of periodic solutions 
implies that the boundary problem for the system (2) with 
boundary conditions 

Zi(t o ) = 0, z(r Q +0) = P(h)z(t +T-0) (4) 
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is solvable for small fj, > 0. Here t is an unknown parame- 
ter, 

/ 1 

P(ji) = -r(jj) 

V £ 2n -2 

where E 2n -2 is the unit matrix of the corresponding size. 
Suppose that the considered boundary problem is solvable 
for fi = 0. Without loss of generality, we assume that the 
corresponding values of t and z(r a ) are zeros. 
Denote 

A= (a tj ) = (A 1 ,...,A 2n ) = (A u A 2 ,A tan ) = 

B = col(6i, . . . , b 2n ) = (5) 
& A m + ^ + 0-0,9 + 0^0, ^)\ Z0=Z ^ ; 
- col(0, /oi, 0, / 2 (0, 0,0),..., 0, /„(0, 0, 0)). 

Now we find the conditions, sufficient for existence of 
a family (t ,z ) = (co1(t Gu), 0, yoo(v), Zo M n{^))) solving 
the boundary problem for the system (2) with the bound- 
ary conditions (4) for small positive /x. It follows from the 
definition of grazing that the inequality 2/oo(m) > must be 
true for positive values of /x. 

Note that b\ — in the considered case. This statement 
follows from the Taylor formula for x\. 

Let ek be the k -th coordinate vector, k = 1, ...,2n, 
Etan = col(e3, . . . , e 2n ) . Due to the implicit function the- 
orem the existence of smooth functions To(^), yoo(^) and 
zo, tan (a*) in the neighborhood of /i = can be provided by 
the condition 



det L y ^ 0; 

L „ = [ (E - A)F a , A 2 + e 2 , A tan - E tan 



(6) 



Here F = F(0, 0,0). The derivatives of r (^), j/oo(m) ancl 
zo,tan{(J') satisfy the system 



L,col(r'(0),^ (0),zU„(0)) 



-B. 



So, to provide the inequality 2/oo(a0 > f° r positive values 
of fi, we must have 



[Ly'Bh < 0. 



(7) 



For the single-degree-of-freedom case the conditions (6) and 
(7) may be simplified. The matrix L y takes the form 

~~ fl 12/oi °12 

(l- 02 2 )/oi a 22 + f/r(0) 



The determinant of this matrix equals to 

A, = -ai 2 /oi(l + l/r(0)). 

So the condition (6) may be rewritten in the form a\2 7^ 0. 
The inverse matrix Ly 1 takes the form 

J_ fa 22 + l/r(0) -aw 
A y \(a 22 - l)/oi -012/01 

and the condition (7) takes the form 

b 2 < 0. 



(8) 



Let us discuss, when the family of periodic solutions, 
corresponding to grazing, persists for after-grazing (nega- 
tive) values of fi. In this case, the boundary problem for 
the system (2) with the conditions 

2/i (r ) = 0, 2(70) = z(t + T), xi(r ) > (9) 

must be solvable. 

Without loss of generality, we may suppose that it is 
the initial conditions t = 0, z(0) = 0, which correspond 
to the solution of the boundary problem (9) for the system 
(2), /i = 0. Then the existence of the smooth functions 
r(/i), xi(fi), Zo t t an (fJ>), giving a solution of the considered 
boundary problem may be provided by the condition 

det L x ^ 0; L x = (A 1 - e u (E - A)F ,A tan - E tan ) 

(10) 

Since foi > 0, the condition (10) may be rewritten in the 
simpler form det(^4 — E) ^ 0. 

To provide the condition x i (fJ.) > for negative values 
of jU it suffices to suppose that 



[L^B]! > 0, 



(11) 



where [L~ 1 B]i is the first element of the vector L~ 1 B. 

In the s.d.f. case the condition (II) may be rewritten in 
a simpler form. If n = 1, 



an - 1 



-012/1 



01 



021 (1 - a 22 )/oi 
The determinant of this matrix equals to 
A, = - det(i4 - E)f 01 



and 



L -\ _ J_ ((A - 022) /01 012/01 
x A x \ -a 2 i an - 1 
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and the condition (11) is equivalent to the following one: 
a 12 b 2 det{A - E) < 0. (12) 

If the conditions (8) and (12) are satisfied simultaneous- 
ly, 012 dct(A - E) > 0. 

Finally, the following two scenarios are possible (the sim- 
ilar results for the square root type mapping are given in 
[11])- 

1. A smooth family ip(t, /x) of periodic solutions exists for 
[i € (/!_,/!+) 9 0. The number of impacts increases 
by 1 as the parameter passes through the bifurcation 
value. This case is called continuous grazing bifurca- 
tion. This bifurcation takes place if (7) and (11) are 
satisfied. For the s.d.f. case this corresponds to the 
case b 2 < 0, ai 2 det(A - E) > 0. 

2. A smooth family <f(t,n) of periodic solutions exists 
for j« G [0,/i + ). These solutions disappear as fi pass- 
es through the bifurcation value. This case is called 
discontinuous grazing bifurcation and corresponds to 
inequalities (7) and 

[L-'Bl < 0. (13) 

For the s.d.f. case the inequality (13) corresponds to 
the case b 2 < 0, ai 2 det(A - E) < 0. 

Remark 4. We mention the following classification of 
periodic solutions in vibro-impact systems, introduced by 
Peterka [1], We say that the periodic solution x(t) of a T 
periodic vibro-impact system is of the (m, n) type (some- 
times the notation m/n is used) if it is of the period nT with 
exactly to impacts over the period. Some bifurcations imply 
the changes in the Peterka pattern. For example, the non- 
degenarate period doubling corresponds to the transition 
(to, n) — > (2m, 2n). In this sense the continuous grazing bi- 
furcation corresponds to the transition (m, n) — > (m± 1, n). 

The bifurcation theory approach to this problem has 
been presented by A. P. Ivanov [48] for the s.d.f. case. The 
following classification has been given: 

1. If Oi2 < (zi2 in the notation of the quoted article) 
then for a positive value of the parameter \i one of the 
multipliers of the given stable periodic orbit turns to 
unity and this orbit collides with the unstable one and 
disappears, i.e. the saddle-node bifurcation scenario 
takes place. 



2. If oi2 > then for a positive value of the parameter 
fi one of the multipliers of the given stable periodic 
orbit turns to —1 and the period doubling bifurcation 
takes place. In this case several periodic solutions of 
different periods may coexist. 

Note that the validity of (8) depends on the choice of 
sign of fi and, if the considered periodic orbit is stable, we 
have det(A — E) > 0. Hence, the existence of a continuous 
grazing family depends on the coefficient 012 only. 

4. Nonhyperbolic chaos in a neighborhood of the 
grazing bifurcation 

Sometimes it may be difficult to apply the Devaney's 
definition of the chaotic invariant set to a multi-dimensional 
VIS. The main trouble is to check the hyperbolicity of the 
obtained invariant set. As we are going to show, one of Lya- 
punov exponents of the periodic solution in the neighbor- 
hood of grazing is large and positive, another one is large 
and negative, and nothing can be said about other ones, 
except they are relatively small with respect to the first 
two. Moreover, these exponents are very sensible to the pa- 
rameters of the system and the choice of the mathematical 
model. 

On the other hand, there are naturally defined stable 
and unstable manifolds of near-grazing periodic solutions 
(both of the dimension 1), corresponding to and the central 
manifold (see the Subsection 4.1 for definition) of the codi- 
mension 2. From the mechanical point of view, this means 
that the normal component x\ of the vector x is "fast" and 
the tangent component x tan is " slow" . 

4-1- Invariant manifolds 

Let U be a domain in the Euclidean space M m . Consider 
a diffcomorphism S G C 1 (U -> V C R m ). Suppose there 
exists an fixed point x* G U of the mapping S. We can 
assume without loss of generality that x* = 0. 

Consider the constants 

< Ai < ,ui < A 2 < fi 2 < A 3 < us, 

fj,! < 1, A 3 > 1 and a decomposition E s © E u © E c = R"\ 
DS(0)(E' J ) = E a for all a G {s, u, c}, such that the follow- 
ing conditions are satisfied. Here [i\ is the degree of contrac- 
tion over the stable space and A3 is the degree of expansion 
over the unstable space. 
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1. Those eigenvalues A of the matrix DS(0), which cor- 
respond to the vectors of the space E s , satisfy the 
estimate Ai < |A| < p\. 

2. Those eigenvalues A, which correspond to the vectors 
of the space E c , satisfy the estimate A2 < |A| < p 2 . 

3. Those eigenvalues A, which correspond to the vectors 
of the space E u , satisfy the estimate A3 < |A| < ^3. 

The spaces E s , E u and E c are called stable, unstable 
and central space respectively. Denote their dimensions by 
m s , m u and m c . Also, we consider the spaces 

E CS =E S ®E C , E CU = E U ®E C . 

Following the notations of [44], we call them central stable 
and central unstable spaces. Let 

m cs = m s + m c , m cu = m u + m c . 

The following result is a corollary of the reduction prin- 
ciple [46] (see also [43] and [47]) and an analogue of the 
well-known Perron theorem for a hyperbolic hxed point. 
Theorem 1. Let S G C 1 ^ -> R m ), G U be a fixed point 
of the mapping S . Then there exist C 1 - smooth embeddings 

b cu : R mC " — >• M. m , such that the following statements hold 
true. 

1. b s (0) = b u (0) = b cs (0) = b cu {0) = 0. 

2. Let 6 CT (K m ") = W? oc , a G {s,u,cs,cu}. These mam- 
folds (we call them local stable, unstable, central stable 
and central unstable respectively) are locally invariant 
in a neighborhood ofO. 

3. The tangent spaces T (VF^ C ) to Wf oc (0) at the origin 
coincide with the spaces E a , a G {s,u,cs,cu}. 

4- For any e > there exists positive numbers a, p and 
5, such that 

\S k (z)\ < (a + e)( Pl +e) k \z\ 
if * e W t s oc , \z\ <5,keN: 
\S k (z)\ < (a + e)(\ 3 -er k \z\ 
¥ zeW^ ocl \z\<5,-keN: 
\S k (z)\ < (a + s)(p 2 +s) k \z\ 
if z G Wf s c , \z\ < S, k G N 
while \S k (z)\ < p; 
\S k (z)\<(a + e)(X 2 -e)- k \z\ 
if z G Wf" c , \z\ < 5, —k G N 
while \S k (z)\ < p. 



5. There exist numbers a > and g > 0, such that for 
any z G B a (0) \ W^ c (0) there is a natural number 
k + such that \S k+ (z)\(pi — e) k+ > g, and for any 
x G Bs(0) \ W^ c (p) there is an integer fc_ < such 
that \S k -{x)\{\ 3 -e) k - > g. 

Remark 5. The stable and the unstable manifolds cor- 
respond to " fast" coordinates and the central manifold cor- 
responds to "slow" coordinate. 

We extend the obtained local manifolds to global ones, 
supposing 

W°(0) = {S k (z) :kEZ,zeW? oc }, 

o G {s, u, cs, cu}. 

The behavior of iterations S k (p) of a point p is the fol- 
lowing: 

1. S k (p) tends to zero exponentially as k — > +00 if 

P eW s (0); 

2. S k (p) tends to zero exponentially as k — > — 00 if 

per(O); 

3. S k (p) does not tend to zero too quickly (this is con- 
trolled by the parameter p 2 ) as k — > +00 if 

p g W cu (0); 

4. S k (p) does not tend to zero too quickly (this is con- 
trolled by the parameter A 2 ) as k — > — 00 if 

p G W cs {0). 

4-2. The separatrix 

Let us consider the system (2), assuming the Condition 
2 is satisfied, i.e. there is a grazing family of periodic solu- 
tions. Fix p, 9 > 0. Denote 

I> - {z° G A : 3h G [-T, T] : Zl (t u -9, z°, p) = 0} 

(Fig. 4). 

Lemma 2. There exists a neighborhood Uo of zero such that 
if the positive parameters p and 9 are small enough, the set 
r M ,e Pi Uq is a surface of the dimension 2n — 1, which is 
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Fig. I 



the graph of the C 2 smooth function x\ — 7 ^fiixtaniV) ■ 
Moreover, 

7 ^ (z)=y K /i(-^,M) + ^ (z) )' (l4) 

where ; y l ^g is a C 2 smooth function such that 7^g(0) = 0. 

Remark 6. This surface T^g plays an important role 
in our later reasonings. We show that corresponds to the 
points of non-smoothness of the Poincare mapping and, at 
the same time, to the solutions with a degenerate impact. 
The invariant manifolds bend, intersecting the surface P^,0- 

Proof. Take a point £ G r^,e- Let the instant to be such 
that zi(to, —0, Q = 0, s = t — to, 

Zi(t + O,-0,C) = col(xi(t),yi(t)). 

Let us show that if to is close enough to —9, we may take 
s o > |*o + 9\ so that the function x\{to + s) does not have 
zeros on [— so,so], except s — 0. Otherwise, there exists 
a sequence £q — s> — (suppose without loss of generali- 
ty, that > —6 and the sequence decreases), a sequence 
t\ G [— 0, £q) and one, consisting of solutions, uniformly 
bounded on the segment [— 0, ij]: 

z*(t) = col(zf(t),...,4(t)) = 

such that (*o) = 0, x\{t\) = (Fig. 5). Also, there exist 
instants if £ (*i)*o)i sucn that ^1(^2) = an d instants 
*a £ (*2i*o) Such that *i(*3) = °- Moreover, i| -> -9, 
xl(tf) -> 0, xf(t§) -> 0. Then 

This contradicts to (3). 




Then for all s € [— so, So] the function x\(to + s) can be 
represented as series 

xi(t + s) =X 2 s 2 +X 3 s 3 + ... (15) 

Differentiating (15), we obtain that 

ii(to + s) = 2X 2 s + 3A 3 s 2 + . . . . 

On the other hand, 

X 2 = x 1 (to + 0)/2->f 1 {-9,z ltt e,li) 
as t -t -0. Then 

xi(-9) = h{-9, z^g, n)(to + 9) 2 {l + o(l))/2; 

yxC— 6>) = /i(-M/», fl ,/*)(to + 0)(i + o(i)). Since yi = 

the formula (14) is true. ■ 

Take a small parameter <; > such that the sets 

F = {z e A : ||z-^,<?|| < ?} c C7 , 
{(x,y) G V : xi < 7 M ,e(a; t(m ,j/)} ) 
{(x,y) G V : xi > Jfj,,e(xtan> y)} 

are correctly defined and nonempty. 

Consider the matrix A, defined by the formula (5). Let 
Ao — det A. Denote the elements of the matrix A by ay 
and ones of the matrix A -1 by ay. Denote the columns of 
matrices A and A 2 by A, and A 2 respectively, the strings 
of the matrix A~ x by Aj. Later on, we assume that 

«X2 > 0, ^axfcafe2 < -ax2- (16) 

As we show later, these conditions allow us to control the 
bent of invariant manifolds, intersecting the surface P^e- 

4-3. The Jacobi matrix 

Note that all the mappings S^ t g, corresponding to the 
same value of /1, are conjugated. Fix a number /1 > 
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and a solution z°(t) = col(ar?(i), y%(t), . . . , x° n (t), y°(t)) of 
the corresponding system with an impact at to- Suppose 
that the corresponding normal velocity Y" i = — 2/i(*o — 0) 
is nonzero. We consider F i as a small parameter. Denote 
z tan(to — 0). Fix a positive number sq and consider 



the mapping 



S (C) = z(*o + «o, to - s , C, A»), 



defined in a neighborhood of the point £° = z°(to — sq). 
Here we assume that the point C° and the parameter sq 
are chosen so that there exists a neighborhood ft 3 £° such 
that any solution z(t) = z(t,t a — So,z_,At) ( z - £ fi) im- 
pacts once over the segment [to — s ,t + Sq\. Denote the 
corresponding instant by t\ = ti(z-) and the normal veloc- 
ity of the impact by Yi = Yi(z-). Let Z tan = z tan (ti) be 
the tangent component of the solution z(t) at the impact 
instant. Take the numbers s± = s±(z-) so that 

to ± s = *i(^-) ± s±(z_) 

for all z_ G f2. The mapping E is smooth in the neigh- 
borhood of the point £o, let us estimate the Jacobi matrix 
£>E(Co). Denote 



z+ 


= z(to - 


fso) 


= z(t - 


f- so,h 


+ 0,0,rY 1 ,Z tan ,fi), 


X+ 


= x(t 


+ s ) 


= x(to 


+ so,h 


+ 0,0,rYi,Z tan ,fi), 


y+ 


= z(to - 


f so) 


= y(to - 


f so,h 


+ 0,0,rY 1 ,Z tan ,/j,), 


X- 


= x(t 


-so) 


= x(t 


- s ,h 


- 0,0, -Y U Z tan ,fl), 


y- 


= y(to - 


-so) 


= y(to - 


- so,h 


-0,0,-Y 1 ,Z tan ,/j 1 ). 



Similarly, we define the numbers 

X\,±, Vl,±, X± > tan, y±,tan- 



Consider the Taylor formula for z±: 

xi,- = Y lS _ + f 1 {t 1 ,0,~Y 1 ,Z tan ,v)s 2 _/2+ 
pi(s-,h, Yi,Z tan ,/j,)s 3 _; 
yi - = -Yi - f 1 (t 1 ,0,-Y 1 ,Z tan ,n)s- + 
p 2 (s-,ti,Yi, Z tan , fj)s 2 _; 

X-.tan = X t an(h) ~ Vtan{t\ ~ 0)s_ + 
ftan(tl,0,-Y 1 ,Z tan , f l)s 2 _/2 + 

p 3 (s-,ti, Yi,Z tan ,/j,)s 3 _; 

ytan- =ytan(h - 0) - ftan(h, 0, -Yi, Z tan , (J,)s. 

p 4 (s-,ti,Yi, Z tan , pt)s 2 _; 

x h+ = rY lS+ + fi(t 1 ,0,rY u Z tan ,ij,)s%/2+ 

p 5 (s + ,tt, Yi, Z tan ,fi)s\; 

= rY x + fi{t 1 ,0,rY 1 ,Z tan ,fj,)s + + 
Pe(s+,ti,Yi, Z tan , m) s +; 

Xtan+ = X t an(h) + ytan(h - 0)s + + 



ftan(h,0, rYi, Ztan, p)s 2 + /2 + 



p 7 (s + ,tl, Yi, Ztan, lA S^A 



Utan-\- ~~ Vtan (ti - 0) + ftan(h,0, rY 1 ,Z tan ,fj,)s++ 

Ps(s+,h,Yi, Ztan, P)s\. 

(17) 

Here all functions, denoted by the letter p with different 
indices, are C 2 smooth with respect to all arguments except 
s±. Denote 

fok+ = fk(to, 0, rY l, Z ,tan, A 4 ), 

fok- = fk{to, 0, — Yoi, ^O.ton, A 4 )- 
It follows from (17) that 
<9z + 



<9(s+, Yi,Z tan ) 




s+=0,Y 1 =Y 01 , Z tan =Z , tan 

(rYoi 
/oi+ r 

y x + o E 2n -2 j 



=0,Y 1 =Y 01 , Z t an=Z ,tan 


-1 

E 2n -2, 



Here E 2n _ 2 is the unit matrix of the corresponding size, 
X+ = col(y 2 (to - 0), / 02 +, • • • , y„(to - 0), f 0n +), 

X- = col(-j/ 2 (*o - 0), -/02+, • • • , -J/n(*o - 0), -/on+)- 

Denote 



w _ dfk(to,0,yi,Zo,tan,p) 

Ik - 



2/1=0 



10 



Clearly, ds + /ds- = — 1. Then, similarly to the results of 
the paper [49], we obtain 



dz + 



B = lim , 

s±^o dz 



-r 
/?2i -r 

v Si ;ta „ E 2n _ 2/ 



Here B M<m = col(/3 3 i, . . . (i 2n l), 

/3 2 i = -(/oi+ + r/oi-)/r 01 = 

-(r + 1)^,(1 + o(y i))/y i, 

/%-!! = 0, 

#y i = (/o,+ - /cy-V^oi = (r + 1)# + O(Foi), 
j = l,...,n. 

Note that dct£ = r 2 . 

It is shown at the Figure 4 how a small rectangular 
domain U of initial data is stretched while the time passes 
a fixed segment corresponding to low-velocity impacts of 
solutions (R transfers to R\). 

4-4- Lyapunov exponents 

In this subsection we study the Lyapunov exponents, 
corresponding to the Poincare mapping at a near-grazing 
periodic point. We show that one of these exponent tend 
to +oo as /i — > 0, another one tends to — oo and all others 
are relatively small, compared with the first two (the ratios 
tend to zero as fi — > 0). 

We estimate the larger and the smaller absolute value of 
eigenvalues of D = DS^g(z^g) and ones of small perturba- 
tions of this matrix. The mapping S^g can be represented 
as the composition S fl _g = S 2 _^_g S\_ fl _g, where 

S2MO =z(T- + 0,0, Cm)- 

The Jacobi matrix A^g = DS 2tlJ/i g(Si il _ li g(z l _ li g)) tends to 
A as [i,0 — > 0. Note, that the elements of these matrices 
satisfy the inequalities similar to (16) provided /1 and are 
sufficiently small. 

The matrix B^g — DS\^^g(z IJ ,^g) is of the form 



lim %± 

s±^o dz- 



B 



as /i and tend to 0. 
Then 

D = A^.gB^g = 

(-(r + l)A 2 <j> (l + o(l))/Y , -rA 2 (l + o(l)), 
A 3 + o(l),...,A 2n + o(l)). 



Consequently, if a\ 2 7^ and if fi is small, one of the eigen- 
values of the matrix D is A + = — (r + l)ai2^o(l + °(l))/^o- 
The corresponding eigenvector u + equals to A 2 + o(l). The 
eigenvalue A + is of the multiplicity 1, the linear space, cor- 
responding to other eigenvalues, tends to the hyperplane 
7Ti, given by the condition x\ — as fi — > 0. Since a\ 2 ^ 0, 
the vector u + is out of 7Ti. Note that 

dct D = dct dct B^g = (r 2 + o(l))A . 

The matrix D" 1 satishes the following asymptotic estimate 



D- L = 



( -rAi+o{\) \ 
(r + lMi(l + o(l))/y 
A3 + o(l) 



V A 2n + o(l) / 

It follows from the form of this matrix, one of the eigenval- 
ues of the matrix D~ x equals to 

AI 1 = (r + 1)^12(1 + o(l))/(r 2 y ). 

The corresponding eigenvector satisfies the asymptotical es- 
timate m_ = e 2 + o(l). 

4-5. Intersection of invariant manifolds 

Let 

A + , A_, A3, . . . , X 2n 
the eigenvalues of the matrix DS l _ lt g(z l _ t ^g). Let 

Ai = min(|A+|, |Al 1 |)/2, A 2 = 2 max (\X k \,\X^\). 

k=3,...,2n 

It is only the first column of D, which is unbounded. Con- 
sequently, Ai > A 2 provided \i is small enough. Let E s 
be the line, corresponding to the eigenvalue A_ , E u be the 
one-dimensional space, corresponding to the eigenvalue A + , 
and E c be the space, corresponding to all other eigenval- 
ues. The spaces E s and E s are linear hulls of the vector 
E{ = e 2 + 0(0) and the vector Ef = A 2 + 0(6) respective- 
ly- 

Due to the Theorem 1, for the fixed point z^ t g there exist 
the stable manifold W s the unstable manifold W u (both 
of the dimension 1), the central stable manifold W cs and 
the central unstable manifold W cu (both of the dimension 
2n — 1). The corresponding spaces E a are tangent to these 
manifolds at the fixed point. 
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Now we show that, provided the assumptions of the The- 
orem 1 are satisfied there exists a non-hyperbolic analogue 
of a homoclinic point for the considered system. 
Lemma 3. There is //+ > 0, so small that, provided 

M 6 (Q,M+) ; 

the manifolds W cs and W u intersect transversally at a point 
P z n,e> i- e - there exists a neighborhood V of the point p, 
such that the connected components of sets W cs |~) V and 
W u P| V , containing p, are smooth manifolds and intersect 
transversally at the point p. 

This intersection for the s.d.f. case is illustrated at the 
figure 6. 




Fig. 6. 



Proof. For small positive values of \i the manifold W u 
intersects transversally the surface r Mi g at the point q u . Let 

= S Plt e{q w )- Denote by L u the arc of the manifold W u , 
connecting the points q u and q\. Denote L" = S llj e(L u ). 
The curve W u is not smooth at the point g". The tangent 
lines to W u at the points of the curve L u are close to the 
line directed along the vector A 2 i.e. the direction of the 
vector A 2 is the limit one to all tangent vectors of W u 
as fi, 6 — > 0. The tangent lines to the arc L" are close to 
the linear hull of the vector A\ for small fi. Due to (16) 
the vectors A 2 and A\ are directed to the different sides 
with respect to the hyperplane tt\ and of the space E cs . 
Consequently, provided \i is small, the neighborhood V can 
be chosen so that the arc L" an d the manifold W cs intersect 
transversally. ■ 

4-6. Straightening 

The main goal of this subsection is to introduce the ap- 
propriate coordinates in a neighborhood of a homoclinic 
point. 



Due to the Theorem 1, there exists a diffeomorphism 
h : U^ : 8 — > f2 of the neighborhood U^g of the point z^fi to 
a domain fi, containing the origin in R 2n and endowed with 
coordinates 

C = coi(c s , C, C) = coi(c s , C u , Ci c , ■ • • , CL-2) 

and such that the following conditions are satisfied (here 
we suppose that h is extended to a diffeomorphism of the 
whole Euclidean space). 

1. h(z lh e)=0. 

2. The Euclidean norms of all columns of the matrix 
Dh{z^$) equal to 1. 

3. MWoc) C {col(C s ,0,0) :( s £l} = OC s , 

h(w^ oc ) c {0,0, C) ■ C eR} = OC, 
KWoc) c {col(C s ,C c ,0) :( s el,C c e M 2 "- 2 }, 
MCc) c {col(0,C c ,0 : C e «-,C £ K 2 "~ 2 }. 

4. Dh{z^e)u + = ex = col(l, 0, . . . , 0). 

5. Dh(z fl ^)u^ = e 2 = col(0, 1, 0, . . . , 0). 

6. The point q — h(p) corresponds to the coordinates 
(C|j0, Cg)i Cq > 0- The segment, linking the point q 
with the origin, is a subset of h(W cs |~| U). 

7. Let be the connected component of h(W u f] U), 
containing the point q. Then 

W)c{coi((;c«,g:c"ei}. 

Let S^.e — h o S^ t $ o h . Consider the point q\ with 
the coordinates (Q, 0, 0). Select the positive numbers e s , e u 
and s c and the neighborhoods 

U = {C-col(C s , C,C C ) : IC1 <eMCI <£",IC C I <£ c }; 

u 1 = {( = co\(cx u ,( c ) ■■ ic s -qi <e s , 

|C"I < e", IC1 < e c } 

(18) 

so that the following conditions are satisfied 

1. U a ,Ui C U. 

2. U f)Ui=(/), i.e. e s < C g s /2- 

3. Uif]h(W s ) i.e. e c > 2|C g c |. 
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4. Let III be the projector to the axis 0( s . For any two 
points 

Ci,2 = col(C 1 ' 3> C u J C c )etf 

the following Lipschitz condition is true: 

M^CiJ-n^te)! < ICf - a 1/2. 

So, we have introduced two disjoint sets: a neighborhood 
Uo of the fixed point and a neighborhood U\ of the homo- 
clinic point. The set U\ is constructed so that it does in- 
tersect the local stable manifold. The reductions \i\S tL fi\u i 
must be uniformly contracting. This may be obtained since 
the distance between the homoclinic point and the fixed 
point tends to zero as the parameter fi tends to the bifur- 
cation value (Fig. 7) . 

4-7. Admissible disks 

Here we consider the set of possible approximations to 
central unstable manifolds and study their iterations. 

We define an admissible disk in one of the domains Ui as 
a subset of Ui, which is the graph of a C 1 - smooth function 
C s = <q{C, CO, defined for |C U | < s u , \( c \ < e c and such that 
max \Dt)(C, C)\ < 1- Let Q l (i = 0, 1) be the sets of disks, 
admissible at Ui endowed with the C 1 metric disti. Denote 

e = QoUCi- 

Any admissible disk intersects transversally with the 
stable manifold. There exists an integer mi such that the 
image S™q(Q) of any admissible disk Q does intersect the 
domain Uo- Moreover, the intersection contains an admis- 
sible curve w i.e. the graph of a C 1 - smooth function 
(C S ,C C ) = x(C"), defined for |C| < e u and such that 



max |x'(C)| < 1, 

|r|<e" 



X(0) = ( X S (0),0) 



i.e. any admissible curve intersects the local stable mani- 
fold. There exist positive numbers /io and £o and the integer 
number m 2 such that if /i < /io, the image S™q(vj) of any 
admissible curve w, satisfying the condition |x(0)| < £o, 
contains an admissible curve w\ C U\. Without loss of gen- 
erality we may suppose m 2 to be so large that for any ad- 
missible disk Q the intersection S™g (Q) f] h(W s ) contains 
a point col(C s ,0,0); |C1 < £o- 

Remark 7. Actually, we are going to present a struc- 
ture, which approximates a foliation of central unstable 
manifolds. In a hyperbolic case we could take a disk, whose 
inclination to the stable manifold is not too small (this is 



what we call admissible disk). Then, due to the A - lemma 
the iterations of this disk provide local approximations to 
the unstable manifold. For non-hyperbolic cases the itera- 
tions of these admissible disks may shrink in some direc- 
tions. This is why we develop a more complex technique, 
using the procedure of extension (see the next subsection). 

Select e\ G (0, e c ] so that for any m > mi the images 
S™g(Q) of any admissible disk Q contain disks, admissible 
at the more narrow domain 

V = {C- \C\ <e s AC\ <e u ,\C\ <£?}■ 

Select the number e\ so that the image S™g(Q) of any disk 

Q, admissible at the set Vq, contains a disk, admissible at 
the domain 

{C: IC S -C 9 1 <e s ,\C\ <eMC c -g < £§} 
and one, admissible at the domain 

{(:\C\<e s ,\C\<e u ,\C\<e c 2 }. 




'0... Fig. 7. 

Denote 5 C — mm(e1, e^), m 3 = mi + m 2 . 

4-8. Extensions 

We are going to show that the initial dynamical system 
engenders one at the set of admissible disks. Unfortunately, 
we cannot consider the iterations S^ e {D) of an admissible 
disk D since these iterations may shrink. So, we need to 
"expand" these iterations on every step. 

Consider the sets 

Vo = {C: |C*| <eMC"l <eMC c | <S c /2} 7 
Vi = {(: |C S - Q\ < e s , \C\ < IC C - C q \ < S c /2}, 

do = {C : IC1 < £ S ,|C"| < eMC c | = S c /2], 
di={(: \C ~Q\< £ s , ICI < £ u , IC C - Q\ = 8 c /2}. 

We can take k so large that the following statements 
hold. 
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1. For all i,j G {0, 1} the intersections 

contain components such that if Q G Qi then the 
intersection S^ e {Q) |"| Uij contains an admissible disk 
Gj(Q) C Vj C Uij. 

2. For all Q 1>2 e Q, i = 0, 1 



disti(Gi(Qi),Gj(Q 2 )) < 2 dist i(Oi-Q2). (19) 



3. The mappings Gj are one-to-one; moreover, the im- 
ages of non-intersecting disks do not intersect any 
more. 

Denote by Qr the sets of the disks, admissible in cor- 
responding domains Vj. Clearly, any disk Q G Q contains a 
subdisk from the set Qr ■ Here we construct an embedding 
J ' : Qr -> Q such that 



Let B = {z G K 2n : |z| < s c }. Define the function 
77 : [-£",£"] x _B by formulae: 



disti(J(Qi), J(Q 2 )) <disti(Qi,Q 2 ) 



(20) 



for all Qi ;2 G Q and J(Q) f| V t = Qf]Vi for all Q G Q~, 
i = 0,l. 

Fix a disk Q G Q~, given by the equation 

Let i G {0,1} be such that Q C XJ{. The intersection 
of the disk Q and the surface dy i is a manifold <9q of the 
dimension 2n — 2. Let us join every point of this manifold 
with a point of the set d\j i by a segment of the type 

{C = col(C s ,C",C c ) :C s =Ci,C u = c 2 , 
C c = tc 3 : r G [t_,t+]} 

so that (ci,C2,t_C3) G 8q, (c\, C2, T-1-C3) G <9j7» or vice versa. 
We obtain a disk in [7^, defined by formulae ( s = 772 (C" j C c )- 
However, this disk is not smooth at the points of the set 

Let us smooth this disk in a standard way. Suppose, 
without loss of generality that i = 0. Consider a C°° - 
smooth function <j> : K — > [0, 1], satisfying the following 
conditions. 

1. (f>{s) = if s < 0, </>(s) = 1 if s > 1. 

2. 0'(s) G (0,2) for all s G (0,1). 



m(CX c ), if|C c l<<572; 

(l-0((2|C c |-<5 c )/ ( 5 c ))r, 1 (C,C c ) + 
0((2|C C |-<5 C )/<5 C )7 ?2 (C,C C ), 

if |C C | g [<57M1; 

r; 2 (C,C c ), if IC C |><5 C - 



Clearly, the graph of this function r\ is an admissible disk. 

Now we fix an admissible disk Q C U [j U\ and apply 
the discussed^ extension-smoothing procedure to the disks 
G {Q) and G\{Q). Denote the obtained disks by G a (Q) 
and G\(Q). The mappings Go and Gi, defined in this way 
are continuous and, due to the formulae (19), both of them 
are uniformly contracting. 

4-9. Symbolic dynamics 

Consider the set E, consisting of infinite one-side se- 
quences a = {ak G {0, 1} : k G Z+}. Let us define the 
metric in the set by a standard way 



d(a,b) =J2 2 fe l afe ~ 



(21) 



k=0 



We identify the periodic subsequences of S, which can be 
obtained by infinite repetition of finite sequences 

(a , . . . , a N ), 

(aj G {0, 1}) with these finite sequences. 

For any periodic sequence a, of this type, generated by 
the set {ao, . . . , oat} we put into correspondence the admis- 
sible disk Q a , which is the unique fixed point of the mapping 

G ao ° G ai O . . . O G aN . 

Then the following statements hold: 

QeS«o, QeG ao (Q ai ), Q e G ao o G ai (Q a2 ), . . . , 
Q G G B0 oG B1 o...oG a „. 1 (Q B)r ). 

(22) 

This means that there exists a constant C > such that 
for any k G N if the first k elements of finite sequences 
a = (ao, ai, . . . , ajvj and 6 = (60, &i, ■ • ■ , &/v 2 ) coincide then 

disti (Q„,Q 6 ) <G2- fe . 



14 



Consequently, for any converging sequence {b k } C E, 
consisting of periodic elements, the corresponding sequence 
Q h k converges at Q. Fix an element 

a = (a ,ai, . . . , a;v, • • •) G E 

and denote a k = (do, . . . , ak)- Clearly, a* — > a. Let 

Q a = limQ afc . 

Consider 

K = {Q a :aeZ}. 

Let H be a parameterizing mapping: H{a) = Q a and cr; be 
the adding of i G {0, 1} to the left side of the element a G E. 
For an arbitrary a G E we denote la = cti(o), 0a = <r n (a). 

Let us prove that the disks, corresponding to different 
elements of the space E, are different. Let a, b G E be such 
that, a^=b. Consider the least integer j, such that aj ^ bj. 
If j = 0, the disks Q a and Qb appertain to different sets £/,. 
Otherwise, due to (22), we have 

Qa G G ao oG ai O...D Ga^iQaj), 

Q b eG ho oG bl o...o G b] _ t (Q bj ) = 
G ao oG ai o...oG aj _ 1 (Q bj ). 

Hence these disks do not intersect. 
Lemma 4. For both i G {0, 1} 

o-ioH = H od. (23) 

or, equivalently, Di a = Gi(D a ) for any a G E, i G {0, 1}. 

Proof. Let i = 0, the case i = 1 is similar. Fix a se- 
quence a — {ak : k G Z + } G E and the corresponding disk 

Qa - Qa„ G ao (Q ai ) Ga ° G ai ( Qa 2 ) R • ■ • 

nG ao oG ai o...oG ajv _ 1 (Q a jn-.- 

Then 

Qoa - Qor\Go(Qa )f)G ° G ao (Q ai )f) - ■ ■ 
f]G oG ao oG ai o...oG ajv _ 1 (Q ajv )n--. = G (Q a ). 

m 

The periodic points of the shift mapping are dense in E, 
and this set is transitive i.e. there is a point 

a* = {a* k ,ke Z+} G E, 

whose shifts are dense [35]. Consequently, the periodic ad- 
missible disks are dense in K and there is a dense sequence 
{Q k : k G Z+} G K such that 

G * {Qk) = Qk-i 



for all k G N. Here all the disks Q k correspond to shifts of 
the sequence a£. 

Note, that the similar set of admissible disks exists for 
any mapping J, satisfying (20). 

4-10. The main result 

Recall some notations: S^g is the Poincare period shift 
mapping for the considered vibro- impact system (2), z^^g 
is the fixed point of this mapping, corresponding to the 
grazing family of periodic solutions, 

dz 

A = lim lim — (T + 9 - 0, + 0, z , fi)U= Zll ,, ■ 
e^o+ oz 

We have proved the following statement. 
Theorem 2. Suppose the vibro-impact system (2) satisfies 
the Condition 2 (existence of the grazing family) and the 
elements of the matrix A satisfy (16). Then there exists 
numbers pi* > 0, 9* > such that for any fj, G (0,/i*) 7 
9 G {0,9*) the domains Ui, defined by (18) and the corre- 
sponding domains Vi, i = 0, 1 are such that the following 
statements are true. 

1. There is an integer k and domains Vq C Uq, V\ C U\, 
such that for any disk D, admissible at Uq or U\, (the 
definition of admissible disk is given at the first para- 
graph of the Subsection 4-1), the set S k g{D) contains 
a disk, admissible at Vq and a disk, admissible at V\ . 

2. The set J of Lipschitz (with the constant, equal to 
1) embeddings J : Q~ — > Q is not empty. Here Q~ 
and Q are the sets of disks, admissible at Vi and Ui 
respectively. 

3. For any J G J there exists a continuous embedding 
H : E — > Q, of the space E consisting of one side 
boolean sequences with the metric (21) to the space Q 
with the G 1 metric. This embedding conjugates map- 
pings Oi ( adding ofO or I to a sequence) and mappings 
d in the sense of formulae (23). 

Remark 8. A non-formal reformulation of this state- 
ment is given at the Section 6. 

4- 11. Structural stability 

Let us discuss the robustness of the constructed set 
with respect to small perturbations of the diffeomorphism 
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S = S^ft or, equivalently, the mapping S — S a ,e- Fix a In this section the following VIS is considered: 

number S > and consider a diffcomorphism S such that 9 

X\ + 2px\ + qx\ — x 2 = 0; x 2 + WI2 = a, 

\S(x) -S(x)\ <S, \DS(x)-DS(x)\<5 while ^(t) > 0; (27) 

the Condition 3 is applied if Xi(t) = 0. 
for any x E Uo\jUi. We suppose the number S be so small ^ ^ ^ 

that for any Q G Q the image S M (Q) contains two admis- x (t) — & + 6sin(cj(t + )) 

sible subdisks of the sets Vi. Also, we suppose that X2 \ > u 2 sin(w^ ^oJJ- 

9 Here 6 and y>o are constants. Then the system (24) may be 

I V £ VQl . . . 



iniS^CijCtan) - ni5(C|,Ctan)| < ^ICi - CI I reduced to the equation 



for all values {Q, (tan) G U. *i + 2 pii + = -5 + &sin(u;(t + <Po))- (28) 

Then, similarly to what we have done above, we may 

construct the set K of admissible disks, with the properties, Provided the conditions (25) are satisfied, the Eq. (28) has 

similar to ones of the set K. the periodic solution 



x\(t) = 

a (q — cj 2 )6sin(cj(t + ip )) 2pbto cos(w(t + <po)) 
uj 2 q Ap 2 uj 2 + (q — to 2 ) 2 4p 2 uj 2 + (q — uj 2 ) 2 



5. Example 

Consider a two degree-of- freedom system, given by equa- 

tl0nS Let ft* (a) = y/Ap 2 uo 2 + {q-uj 2 ) 2 a/{quj 2 ). If b = 6* (a), the 

x 1 + 2p±i + gsi - x 2 = 0; x 2 + w 2 .x 2 = a. (24) function x\ (t) is non-negative and there is exactly one im- 
pact instant of the corresponding solution over the period 
Suppose that [0: 27t/cj), which is grazing. The general solution of the Eq. 



q-p 2 =uj 2 >0, a > 0, u> > 0, p 2 + (q-uj) 2 > (25) 



(28) is of the form 



Xl (t) = Cexp(-pt) sm(uj (t + 9)) + 



and there exists k e N, such that , ■>%,.,/ „ „ , , w 9 

(q — uj )bsm(uj(t + (fo)) — 2pbuj cos{uj{t + <p )) 

luo/uj G (fc + l/4,fc + l/2). (26) 4p 2 w 2 + (q - w 2 ) 2 

Denote T = 2^/lo. Suppose that the system (24) is dchncd Hcrc C and 9 are constants. All the (1,1) periodic solu- 
for Xl > 0. The impact conditions are the following. tions of the s y stem ( 27 ) are ones of the boundary problem, 

Condition 3. depending on the parameter t G [0, T) 

1. Let re (0,1]. If x 1 (t )=0, then ( x(t ) = x(t + T) = 0; x(t ) = -rx(t + T); 

{ x(t)>0 forte (t , to +T). ^ j 

a;i(i o + 0)=0, x 2 (t + 0) = x 2 (t - 0), f 

Fix 6 and </? - The Eq. (29) may be considered as a 

ii(t + 0) = -rxi(t - 0), x 2 (t + 0) = x 2 (t - 0). system on the parameters C, and t . Denote 



2. Let ffce instant t be such that Xl (t ) = ±i(t - 0) = *i = Kcsm(2pw/ ^pV + fe-w) 2 ), 

and x 2 (t ) < 0. Consider the solution £(t) of the sec- tq = t + (p - ti, t = t + Lp Q - ti, 

ond equation (24) with initial conditions A = a/ {u 2 q), B = b/ ^4p 2 w 2 + (q - ui) 2 ), 



£(to) = x 2 (to - 0), £(io) = i 2 (t - 0). 



d = Cexp(pt ), tf = + t o . 
The general solution of the Eq. (28) is of the form 



Let I — [to,ti] be the maximal segment such that 

x 2 (t) < for all t e I. Then Xl (t)\i = 0. x ^ = C ^ exp(-p(r - r )) sin( Wo (r - r + 0))+ 

A + J Bsin(wr). 
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The first two equations of (29) may be rewritten as follows: 

C\ sin(wo^) + A + B sin(wT ) = 0; 
Ci sin(wo^) = Ci cxp(-pT) sin(w (?? + T)); 
< Ci(ujo cos(wo^) — _psin(wo^)) + Bujcos(ojt ) = (30) 
-rCi cxp(-pT)(w cos(w (tf + T))- 
psin(cj (^ + T))) - rBbj cos(w(t + T)). 

The second of the Eq. (30) gives two cases: either C\ = 0, 
or 

./ <n exp(pT) - cos(w T) 

COt UdnU = ; — . (31) 

V ; sin(u;oT) 1 ; 

If C\ = 0, then b = b*(a) and function X\{t) has exactly 
one zero (of the multiplicity 2) over the period. Otherwise 
the condition (31) defines uniquely the number 

•d e (0,7t/cj ). 

It follows from the third equation (30) that 

d = B{1 + r)u cos(ut q )/H. (32) 

Here 

H = rexp(-pT)(psm(oj (ti + T))- 
loq cos(wo($ + T))) + psin(wo'i9) — wo cos(cjo??) • 

Substituting the expression (32) to the first equation (30), 
we obtain 

B(sm(ujT ) + D cos(wto)) + A = 0, (33) 

where D = (1 + r)w sm(uio'&)/H. The number To may be 
found from the Eq. (33) if and only if 

B 2 (l + D 2 )>A 2 . (34) 

This inequality implies \b\ > 6* (a). If (34) is false, the sys- 
tem (29) is unsolvable. For b = ±b*(a) the graph of the 
function x\ (t) is tangent to Ox and does not have any oth- 
er zeros over the period. If \b\ > b* (a), then for any fixed 
value of the pair (b, ipo) there are two solutions of the Eq. 
(33) and, respectively, two periodic solutions of the VIS 
(27). These solutions continuously depend on b and tpo and 
have a single impact over the period. The corresponding 
velocities tend to as b — > ±b*(a). The system 

y = -qx - 2py 



is the homogeneous part for (28). Let $(t) be the funda- 
mental matrix of the system (35), turning to the unity at 
t = 0. The trace of $(T) equals to 2 cos(aj T) exp(-pT), 
and the determinant equals to cxp(— 2pT). Due to the in- 
equality (26), the conditions of the Theorem 2 are satisfied 
and there exists a set of admissible disks, corresponding to 
the VIS (27), which can be described by means of symbolic 
dynamics. 

6. Discussion. 

Here we reformulate the result of the Theorem 2 on a 
less formal level. We try to find a "visible" analogue of the 
chaotic dynamics in a neighborhood of an invisible chaotic 
invariant set. First of of all we construct a non-hyperbolic 
homoclinic point. 

For this we need conditions (16). These conditions have 
a clear geometric sense. The inequality ayi > implies that 
the unstable manifold W u exists and intersects the grazing 
surface T^j provided fi and 9 are positive and sufficient- 
ly small. The unstable manifold bents at the point of the 
intersection. If 

2n 

^ aik0.k2 < —0,12 
J = l 

the manifold W u turns so that it intersects W cs transver- 
sally. 

This does not imply the existence of a Smale horseshoe. 
What we can expect is to find a symbolic dynamics on leaves 
of a foliation over a Cantor set. However, even this is not 
always possible. 

A central stable manifold is not as "good" as an ordinary 
stable manifold. For example, it may be non-unique even for 
a fixed point [44] . Typically, a central stable bundle over an 
invariant set is not integrable [43] i.e. we can not find a 
corresponding continuous foliation, like what we can do in 
the hyperbolic case. Moreover, the iterations of a central 
stable (or a central unstable) manifold may shrink both in 
positive and negative direction. So, we have to be able to 
expand iterated leaf on every step. This is why we need the 
mapping J from the statement of the theorem. 

If we can reconstruct the whole central stable (or central 
unstable) leave by any small plaque of this disk (this is the 
so-called dynamical coherence [45] ) , then we get a naturally 
defined expanding mapping J and, due to the Theorem 2, 
a symbolic dynamics on the set of central unstable leaves. 
For example, this happens, if the considered central unsta- 
ble manifolds are de facto ordinary unstable manifolds (a 
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hyperbolic case) or, more generally, there is a Lyapunov 
function, which does not allow the central unstable leaves 
to shrink. For the last case, we do not need any expanding 
mappings at all. 

In the general case, we offer a sample of expanding map- 
ping (Subsection 4.8) drawing straight lines through the 
boundary points of the "narrow" disk. The expansion pro- 
cedure allows us to construct a foliation. So, iterating and 
expanding disks on every step, we obtain a complex dynam- 
ics on the set of disks, close to the central unstable bundle. 
We have constructed a model of a slow-fast system with a 
chaotic dynamics for the fast variable and unknown behav- 
ior of the slow one. 

Compare the Devaney's definition of chaos with one of- 
fered in our paper. 

Remark 9. 

1. The Devaney's definition of chaos requests hyperbol- 
icity. The introduced model does not. 

2. The Devaney's definition describes the dynamic of 
point. Our model describes the dynamic of approx- 
imations to central unstable leaves. 

3. Both the considered models are structurally stable. 

7. Conclusion. 

An impact oscillator, described by the Newtonian mod- 
el of impacts has been studied. There exist many results 
on grazing bifurcations for such systems (see [10-12] for re- 
view). It has been shown that different types of chaos can be 
observed for near-grazing values of parameters. The explicit 
conditions, providing existence of a complex dynamics have 
been given. The geometry of the obtained strange attractors 
has been described. 

However, especially in the case of several degrees-of- 
freedom, the structure of these attractors may be very sensi- 
tive to the parameters of the system. This happens because 
the ratios of the Lyapunov exponents of this solution are 
big in a non-degenerate case. 

In order to describe the robust properties of attractors, 
we offer a new approach treating the near-grazing dynamics 
as non-hyperbolic. We use a technique of central unstable 
manifolds. A dynamical system on a set of embeddings of 
the local central unstable manifold of a periodic point is in- 
troduced. For this system we find an invariant set described 
by the symbolical dynamics. 



The conditions of existence of this set (Theorem 2) do 
not require hyperbolicity of the periodic motion. So, the 
obtained structure is less sensitive to the changes of pa- 
rameters when the strange attractor, which may disappear 
if one of the Lyapunov exponents becomes zero. Roughly 
speaking, some traces of the strange attractor persist, even 
when the attractor itself has already disappeared. 

From the mathematical point of view, we do not use any 
common assumptions from the theory of partial hyperbol- 
icity, like integrability of the central unstable bundle, etc. 
There may be no continuous central unstable foliation in 
the considered case. What we have constructed, is not a 
classical "dynamics of leaves" [43]. However, all the "ad- 
missible disks" are close to the central unstable bundle and 
can be considered as approximations to local central unsta- 
ble manifolds. So, the result of the Theorem 2 describes a 
new type of chaotic dynamics. 

Grace to the work by Dankowitz and Nordmark [7] , the 
results of the current paper may be transferred to the "soft" 
models of impacts. 

The obtained model of near-grazing dynamics is struc- 
turally stable and can be observed in experiments. 
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